Four small puzzles that Rosetta doesn't solve

A complete macromolecule modeling package must be able to solve the simplest structure prediction problems. Despite recent successes in high resolution structure modeling and design, the Rosetta software suite fares poorly on deceptively small protein and RNA puzzles, some as small as four residues. To illustrate these problems, this manuscript presents extensive Rosetta results for four well-defined test cases: the 20-residue mini-protein Trp cage, an even smaller disulfide-stabilized conotoxin, the reactive loop of a serine protease inhibitor, and a UUCG RNA tetraloop. In contrast to previous Rosetta studies, several lines of evidence indicate that conformational sampling is not the major bottleneck in modeling these small systems. Instead, approximations and omissions in the Rosetta all-atom energy function currently preclude discriminating experimentally observed conformations from de novo models at atomic resolution. These molecular "puzzles" should serve as useful model systems for developers wishing to make foundational improvements to this powerful modeling suite.Comment: Published in PLoS One as a manuscript for the RosettaCon 2010 Special Collectio

Increasing diterpene yield with a modular metabolic engineering system in E. coli: comparison of MEV and MEP isoprenoid precursor pathway engineering

Engineering biosynthetic pathways in heterologous microbial host organisms offers an elegant approach to pathway elucidation via the incorporation of putative biosynthetic enzymes and characterization of resulting novel metabolites. Our previous work in Escherichia coli demonstrated the feasibility of a facile modular approach to engineering the production of labdane-related diterpene (20 carbon) natural products. However, yield was limited (<0.1 mg/L), presumably due to reliance on endogenous production of the isoprenoid precursors dimethylallyl diphosphate and isopentenyl diphosphate. Here, we report incorporation of either a heterologous mevalonate pathway (MEV) or enhancement of the endogenous methyl erythritol phosphate pathway (MEP) with our modular metabolic engineering system. With MEP pathway enhancement, it was found that pyruvate supplementation of rich media and simultaneous overexpression of three genes (idi, dxs, and dxr) resulted in the greatest increase in diterpene yield, indicating distributed metabolic control within this pathway. Incorporation of a heterologous MEV pathway in bioreactor grown cultures resulted in significantly higher yields than MEP pathway enhancement. We have established suitable growth conditions for diterpene production levels ranging from 10 to >100 mg/L of E. coli culture. These amounts are sufficient for nuclear magnetic resonance analyses, enabling characterization of enzymatic products and hence, pathway elucidation. Furthermore, these results represent an up to >1,000-fold improvement in diterpene production from our facile, modular platform, with MEP pathway enhancement offering a cost effective alternative with reasonable yield. Finally, we reiterate here that this modular approach is expandable and should be easily adaptable to the production of any terpenoid natural product

Expression of Fraser syndrome genes in normal and polycystic murine kidneys

BACKGROUND: Fraser syndrome (FS) features renal agenesis and cystic kidneys. Mutations of FRAS1 (Fraser syndrome 1)and FREM2 (FRAS1-related extracellular matrix protein 2)cause FS. They code for basement membrane proteins expressed in metanephric epithelia where they mediate epithelial/mesenchymal signalling. Little is known about whether and where these molecules are expressed in more mature kidneys. METHODS: In healthy and congenital polycystic kidney (cpk)mouse kidneys we sought Frem2 expression using a LacZ reporter gene and quantified Fras family transcripts. Fras1 immunohistochemistry was undertaken in cystic kidneys from cpk mice and PCK (Pkhd1 mutant) rats (models of autosomal recessive polycystic kidney disease) and in wildtype metanephroi rendered cystic by dexamethasone. RESULTS: Nascent nephrons transiently expressed Frem2 in both tubule and podocyte epithelia. Maturing and adult collecting ducts also expressed Frem2. Frem2 was expressed in cpk cystic epithelia although Frem2 haploinsufficiency did not significantly modify cystogenesis in vivo. Fras1 transcripts were significantly upregulated, and Frem3 downregulated, in polycystic kidneys versus the non-cystic kidneys of littermates. Fras1 was immunodetected in cpk, PCK and dexamethasone-induced cystepithelia. CONCLUSIONS: These descriptive results are consistent with the hypothesis that Fras family molecules play diverse roles in kidney epithelia. In future, this should be tested by conditional deletion of FS genes in nephron segments and collecting ducts

Enhancing Production of Bio-Isoprene Using Hybrid MVA Pathway and Isoprene Synthase in E. coli

The depleting petroleum reserve, increasingly severe energy crisis, and global climate change are reigniting enthusiasm for seeking sustainable technologies to replace petroleum as a source of fuel and chemicals. In this paper, the efficiency of the MVA pathway on isoprene production has been improved as follows: firstly, in order to increase MVA production, the source of the “upper pathway” which contains HMG-CoA synthase, acetyl-CoA acetyltransferase and HMG-CoA reductase to covert acetyl-CoA into MVA has been changed from Saccharomyces cerevisiae to Enterococcus faecalis; secondly, to further enhance the production of MVA and isoprene, a alanine 110 of the mvaS gene has been mutated to a glycine. The final genetic strain YJM25 containing the optimized MVA pathway and isoprene synthase from Populus alba can accumulate isoprene up to 6.3 g/L after 40 h of fed-batch cultivation

A Kinetic Model of Trp-Cage Folding from Multiple Biased Molecular Dynamics Simulations

Trp-cage is a designed 20-residue polypeptide that, in spite of its size, shares several features with larger globular proteins. Although the system has been intensively investigated experimentally and theoretically, its folding mechanism is not yet fully understood. Indeed, some experiments suggest a two-state behavior, while others point to the presence of intermediates. In this work we show that the results of a bias-exchange metadynamics simulation can be used for constructing a detailed thermodynamic and kinetic model of the system. The model, although constructed from a biased simulation, has a quality similar to those extracted from the analysis of long unbiased molecular dynamics trajectories. This is demonstrated by a careful benchmark of the approach on a smaller system, the solvated Ace-Ala3-Nme peptide. For the Trp-cage folding, the model predicts that the relaxation time of 3100 ns observed experimentally is due to the presence of a compact molten globule-like conformation. This state has an occupancy of only 3% at 300 K, but acts as a kinetic trap. Instead, non-compact structures relax to the folded state on the sub-microsecond timescale. The model also predicts the presence of a state at of 4.4 Å from the NMR structure in which the Trp strongly interacts with Pro12. This state can explain the abnormal temperature dependence of the and chemical shifts. The structures of the two most stable misfolded intermediates are in agreement with NMR experiments on the unfolded protein. Our work shows that, using biased molecular dynamics trajectories, it is possible to construct a model describing in detail the Trp-cage folding kinetics and thermodynamics in agreement with experimental data